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1. Introduction 

The spacelike pion electromagnetic form factor F n (t) is of interest for studying the onset of 
perturbative QCD (pQCD). At leading order (LO), the expression for F K (t), calculated in pQCD is 
given by 



pert 



Q 2 



t(=-Q 2 )<0 



(l.l) 



while the next-to-leading order (NLO) correction with the MS-renormalization scheme and asymp- 
totic DAs reads [Si 
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where f n = 130.4 MeV is the pion decay constant and a s (fL 2 ) the strong coupling at the renor- 
malization scale /I 2 . The quantity j8o = 11 — 2«y /3 denotes the first coefficient in the perturbative 
expansion of the j8 -function, rif being the number of active flavors. 

There is an interplay of perturbative with soft, nonperturbative processes especially in the 
intermediate Q 2 region as a result of which the asymptotic regime sets in quite slowly in the case 
of the pion form factor. Therefore, it remains an open question as to at what value of Q 2 do 
the nonperturbative contributions become negligible so that the perturbative QCD description of 
the form factor becomes reliable. Several nonperturbative approaches have been proposed for 
the study of the spacelike pion form factor [Q, [7| [?[ 10, 11, 12|. On the experimental side, 
measurements of the spacelike pion form factor at various energies are available, the most recent 
data coming from the J LAB, Ji3| , 14]. Further, we now have more accurate information on the 



phase [15] and modulus [16] of the pion form factor on the unitarity cut. In this paper, we perform 
an analytic continuation from the timelike to the spacelike region using in a most conservative way 
the available information on the phase and modulus on the unitarity cut, and also the spacelike 



information available. Using a mathematical formalism discussed in []17j, |18Q , we find stringent 
upper and lower bounds at different values of spacelike momenta, which allows us to find a lower 



limit for the onset of the QCD perturbative behavior [19 ]. 



2. Formalism 

Our formalism requires the knowledge of the phase below t m and an integral over the modulus 
squared from t m to oo. We relate the phase of pion the form factor with that of the associated %% 
scattering amplitude via the Fermi Watson theorem. In this case, we consider the relation 

Arg[F(f + is)) = 81 (t), 4M 2 K <t< t in , (2.1) 

where 5f (?) is the phase shift of the P wave of nn elastic scattering and t m = [M n + M a ) 2 is the 
upper limit of the elastic region, which corresponds to the first important inelastic threshold due to 
the (on pair. We use the recent experimental data on the modulus up to \ft = 3GeV [ 16]. Above 
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this energy, we make conservative assumptions and obtain a rather accurate estimate of an integral 
of modulus squared from t m to infinity. More precisely, we assume the following condition, 



1 

nJt 



dtp{t)\F n {t)\ l =I, (2.2) 



where p{t) is a suitable positive-definite weight, for which the integral converges, and the number 
/ can be estimated with sufficient precision. The optimal procedure is to vary p(t) over a suitable 
admissible class and take the best result. In principle, a large class of positive weights, leading to 
a convergent integral for compatible with the asymptotic behavior ( |1.1[ ) of the pion form 

factor, can be adopted. In our calculations, we consider an expression of the form 

Pa(t) = L 0<a<2 (2.3) 
t 

We use additional information inside the analyticity domain namely the normalization Fj^O) = 
1 and the pion charge radius, F^(0) = with (r 2 ) varied within reasonable limits [H]], 



and the values of the form factor at some spacelike values F(t n ) where t n < Q13| , |14Q . In this 
paper, we derive rigorous upper and lower bounds on F K (t) in the region t < 0, for functions F n (t) 
belonging to the class of real analytic functions in the ? -plane cut for t > 4M 2 defined by all the 
inputs specified. 

For solving the problem, we apply a standard mathematical method discussed in detail in 



[17, 18]. We transform our problem via a conformal map, cast the integral equation into a canonical 



form and derive a determinant (see ref. [ 19] for more details) which is central to our investigations 



for obtaining bounds on F n (t) in the spacelike region. 



3. Inputs 

In the elastic region t < t m we use the phase shift parameterization determined recently with 



high precision from Roy equations applied to the kk elastic amplitude given in | |15| ] (see ref.[]19|] 
for more details). Above t m we choose 8 as a continuous function, sufficiently smooth which ap- 
proaches asymptotically %. The results are independent of the choice of 8(t) above t m , as discussed 



in detail in [|18fl. For the calculation of the integral defined in (|2.2|), we use the BaBar data [ |16| ] 
from = 0.917 GeV up to y/t = 3 GeV, and have taken a constant value for the modulus in the 
range 3 GeV < y/t < 20 GeV, continued with a l/t decrease above 20 GeV. It may be noted that 
our estimates are based on very conservative assumptions of the input quantities which makes our 



procedure very robust. In our analysis, we consider the weights of the form given in ( |2.3| ). The 
values of / corresponding to several choices of the parameter a are given in Table 1 of ref. [|19|], 
where the uncertainties are due to the BaBar experimental errors. We find from our analysis that 
the best results come from an optimal choice of the weight corresponding to Pi/2(?)- We use also 
as input, 

(4) =0.43 ±0.01 fm 2 , F(-2.45GeV 2 ) =0.167 ±0.010t{j;g$, (3.1) 

respectively for the pion charge radius [^, 21 ] and the spacelike datum [|13|, 14]. 



3 



Bounds on the spacelike pion electromagnetic form factor 



I. Sentitemsu Imsong 



4. Results 

The main results emerging from our analysis are presented in Figs. [T] and ^[ In the figures, 
the white band conesponds to the bound obtained by using only the central values of the inputs 
while the grey bands are obtained from the errors associated with the inputs. The error bands have 
been obtained by adding quadratically the errors produced by the variation of the phase, the charge 
radius, the integral / for a = 1/2 from Table 1 of ref. [|19|], and the spacelike datum. We find that the 
greatest contribution to the size of the grey domain is the experimental uncertainty associated with 
the spacelike value (3T). In Fig. [j], we compare our constraints with some of the data available 



from experiments (see [19| for references). We find that at lower values of Q most of the data 
are consistent with the narrow allowed band predicted by our analysis. In Fig. ||, we compare our 
allowed domain with the pQCD predictions both at LO and NLO and with various nonperturbative 
models. The perturbative prediction to NLO is sensitive to the choice of the renormalization scale, 
and also of the factorization scale in the case when pion DAs different from the asymptotic ones are 
used in the calculation ||5|, |22|] . Several prescriptions for scale setting have been adopted, but there 
is no general consensus on the issue. For illustration, in Fig. |2] we show the sum of the LO and 



NLO terms ( |1.1[ ) and (1.2), obtained with the scale }i 2 = Q 2 and the one loop coupling. This curve 



is compatible with our bounds enlarged by errors only for Q 2 > 6 GeV 2 . We show also several 
nonperturbative models proposed in the literature for the spacelike form factor at intermediate 
region [||, 0, |], ||, 10, 12]. In Ref. the authors applied light-cone QCD sum rules and 
parametrized with a simple expression the nonperturbative correction, to be added to the LO+NLO 
perturbative prediction in the region 1< Q 2 < 15 GeV 2 . We show the sum of the soft correction and 
the perturbative QCD prediction to NLO, evaluated at a scale \i 2 = 0.5 Q 2 +M 2 with M 2 = 1 GeV 2 
as argued in [|]]. The model is quite compatible with our bounds, the corresponding curve being 
inside the small white inner domain for Q 2 > 6 GeV 2 . The model based on local duality [[7|] is also 
consistent with the allowed domain derived here for Q 2 > lGeV 2 . We mention that this model, 



proposed in Q23Q, was recently developed by several authors [|1 1|]. The other models shown in 
Fig. || are consistent with the bounds derived by us at low Q 2 , but are at the upper limit of the 
allowed domain at higher Q 2 . The agreement is somewhat better for the model discussed in [B], 
which is a LO+NLO perturbative calculation using nonasymptotic pion DAs evolved to NLO, with 
a modification of the QCD coupling by the so-called analytic perturbation theory. The AdS/QCD 
model considered in [§] is in fact a simple dipole interpolation, which is valid at low energies but 
seems to overestimate the form factor at larger momenta. The same remark holds for the models 



discussed in [10] and [12], based on QCD sum rules with nonlocal condensates, and the chiral limit 



of the hard-wall AdS/QCD approach, respectively. 



5. Conclusion 

We have derived upper and lower bounds on the pion electromagnetic form factor along the 
spacelike axis, by exploiting in a conservative way the precise information on the phase and mod- 
ulus of the timelike axis as well as the available spacelike data. We have used the method of 
analytic continuation to obtain information on the spacelike region of the pion form factor. Using 
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Figure 1: Allowed domain obtained with the weight Pi/2(f) compared with several sets of experimental 
data. 
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Figure 2: Comparison of the bounds for the weight P\/2(t), with perturbative QCD and several nonpertur- 
bative models. 



the weight P\/2{t) which is the optimal choice as discussed in Sec. ^ we have obtained upper and 
lower bounds enlarged by errors associated with the various inputs entering our analysis. From Fig. 
||[ we can conclude that perturbative QCD to LO is excluded for Q 2 < 7 GeV 2 , and perturbative 
QCD to NLO is excluded for Q 2 < 6 GeV 2 , respectively. If we restrict to the inner white allowed 
domain obtained with the central values of the input, the exclusion regions become Q 2 < 9 GeV 2 
and Q 2 < 8 GeV 2 , respectively. Among the theoretical models, the light-cone QCD sum rules [||] 
and the local quark-hadron duality model [J7]] are consistent with the allowed domain derived here 
for a large energy interval, while the remaining models are consistent with the bounds at low ener- 
gies, but seem to predict too high values at higher Q 2 . To increase the strength of the predictions, a 
reduction of the grey bands produced by the uncertainties of the input is desirable. As mentioned 
in SecM the biggest contribution to the error band comes from the experimental errors associated 
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with the spacelike datum (3.1). As such, in order to increase the predictive power of our formal- 
ism, more accurate data at a few spacelike points, particularly at larger values of Q 2 , would be very 
useful. 
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